Centrality Scaling of the pr Distribution of Pions 



Rudolph C. Hwa 1 and C. B. Yang 1 ' 2 
1 Institute of Theoretical Science and Department of Physics 
University of Oregon, Eugene, OR 97403-5203, USA 
2 Institute of Particle Physics, Hua-Zhong Normal University, Wuhan 430079, China 

From the preliminary data of PHENIX on the centrality dependence of the ty° spectrum in pr at 
midrapidity in heavy-ion collisions, we show that a scaling behavior exists that is independent of 
the centrality. It is then shown that (pr) degrades with increasing 7V par t exponentially with a decay 
constant that can be quantified. A scaling distribution in terms of an intuitive scaling variable is 
derived that is analogous to the KNO scaling. No theoretical models are used in any part of this 
phenomenological analysis. 



PACS numbers: 25.75.Dv 
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In a recent paper we reported on the finding of a 
scaling property of the pr distribution of pions produced 
in heavy-ion collisions that is independent of the colli- 
sion energy. Here we present an extension of that scaling 
property to include centrality variations and show that 
a KNO-type scaling behavior exists over the entire 
range of pr measured. The investigation is primarily 
a phenomenological analysis with no assumptions about 
the hard and soft collisions, nor about the parton energy 
losses. 

Recently, a scaling behavior of the transverse-mass 
spectrum has been reported in . That work was moti- 
vated by color glass condensate and the saturation of the 
gluon density in nuclear collisions. Our investigation has 
no theoretical motivation other than the search for the 
simplest form that can represent the data. The dynam- 
ical origin of the px distribution is complicated. At low 
Pt the statistical model seems to work well, as does the 
hydrodynamical description up to pr =3 GeV/c |^|. At 
high pt hard parton scattering will create jets, which can 
lose energy due to multiple scatterings of partons in the 
dense medium Q . A universal description of the hadron 
distribution over all px is nonexistent, if not meaning- 
less from the point of view of the sectarian nature of the 
dynamical theories that claim validities in different do- 
mains. However, if a universal scaling behavior can be 
found phenomenologically, it can serve as a common goal 
for different dynamical approaches to aim at. 

From the preliminary PHENIX data of ir° pro- 
duced in Au+Au collisions at the relativistic heavy- 
ion collider (RHIC) @, we have the pion distribution, 
(2-Kp T y 1 dN 7r /dp T , at midrapidity for ^/s = 200 GeV 
and for a wide range of centrality that has 9 bins from 
0-10% to 80-92%. To unify the 9 distributions, it is 
necessary to define a scaling variable z. First, we use 
the number of participants, AT par t, to quantify central- 
ity; those numbers for different bins are taken from 0, 
which agree well with those given by PHENIX 0. Next, 
we define, for fixed y/s (at 200 GeV), 



z = PT /K{N), 



(1) 
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FIG. 1: Scaling distribution 4>(.z) showing the coalescence of 
5 centrality bins of the preliminary data from PHENIX on 
7T° production in Au+Au collisions at yfs = 200 GeV [g. 
The points labeled by x are for 7r + . The solid line is a fit 
parametrized by Eq. J^J. 



bin we vary K by plotting the data of (2ttpt)~ 1 dN n / dpr 
in terms of z and adjusting the normalization so that all 
data points lie on a universal curve. That is, we define 



$0) = A(N) K 2 (N) 



1 dN« 
2ttpt dpT ' 



(2) 



where K depends on AT part , for which we use the abbrevi- 
ated notation N = iVp art hereafter. For every centrality 



and find A(N) and K(N) such that $>(z) has no explicit 
dependence on N. That turns out to be possible, as 
evidenced by Fig. 1. For clarity we show only 5 bins of 
centrality in that figure. It is a remarkable property of 
the centrality dependence of the pion spectra that such 
a universal scaling distribution exists. 

The values of K(N) used to obtain the scaling behavior 
are shown in Fig. 2(a) in units of GeV/c. The dependence 
of K(N) on N can be well fitted by 



K(N) = 1.226 - 6.36xlO~ 4 Af, 



(3) 
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FIG. 2: (a) Scale factor K(N) in units of GeV/c. Solid line is 
a fit by Eq. (b) Power- law behavior of the normalization 
factor A(N C ). Solid line is a fit by Eq. (gj. 

such that AT(iV m ax) = 1 at N = N max = 350. The effects 
of the degradation of parton momenta are hidden in this 
formula. Any change of the overall scale of K(N) is triv- 
ial and does not affect the scaling behavior that we have 
found. Although the normalization factor A(N) does not 
have a simple dependence on A, it turns out to depend 
simply on the number of binary collisions N c . The values 
of A(N C ) needed to achieve the scaling $(z) are shown 
in Fig. 2(b) in a log-log plot. They can be fitted by 

A{N C ) = 530A,T - 9 . (4) 

From the tables listed in Refs. 0, , N c and N can be 
related by N c = 0.44./V 1 - 33 . Note that the normalization 
of "3?(z) is set by the most central bin by choosing A(N) = 
1 at N = A max . If A(N C ) were to behave as N' 1 , it 
would suggest that the average multiplicity of pions at 
midrapidity is proportional to N c , which is a variable 
that measures the number of hard collisions. Thus the 
factor N~ - 9 in Eq.(gJ is an indication that the centrality 
dependence of the midrapidity multiplicity scales as N® 9 
from the pp collisions, revealing the effect of suppression 
of pt in the nuclear medium. 

To fit the scaling curve the ir° data are insufficient to 
give us guidance in the small z region, since they do not 
extend below pr=l GeV/c. For < pr < 1 GeV/c, 
we use the tt+ data of PHENIX for 0-5% centrality @ 
shown in Fig. 1. The combined tt° and ir + data can be 
well fitted by 

*(«) = 1200(z 2 + 2)- 4 - 8 (l + 25e- 4 - 52 ), (5) 

which is shown by the solid line in Fig. 1. We can check 
its normalization by evaluating the integral 

f W , , A(N) dN^o 

1= dzz$(z =46.2= r — • (6) 
Jo 2ir d V V ; 



For N — 200, say, this gives dN^a /drj = 149, which com- 
pares satisfactorily to dN c h/dr]/(0.5N) — 3.2 at the same 
N . Since the 7r ± data do not extend into the pr > 2 
GeV/c region, we do not consider them for centrality 
analysis here. 

The exponential term in Eq. (JjjJ is mainly to fit the low- 
z data that contain thermodynamical effects. At high z, 
&(z) behaves as a power law that represents the effects 
of hard collisions and jet quenching. For all z, &(z) is a 
succinct summary of all dynamical effects for all central- 
ities. 

In terms of $(z) it is now possible to have an analytic 
expression of the inclusive distribution of the pions in 
Pt at midrapidity. For convenience, we shall write it in 
terms of the momentum fraction x 

x = pt/K , (7) 

where K is a fixed scale, beyond which no physics of 
interest need be of concern here. We set Kq = 10 GeV/c 
for now, although increasing it later, if necessary, is a 
simple matter. In view of Eq. we thus have 

z = xA(N), A(N)=K /K(N). (8) 

Converting (2irz)~ 1 dN 7T /dz to the x variable, we define 
the corresponding pion distribution to be 

H{ Xl N) = A- 1 {N)<^{x,N), (9) 

where A(N) = A(N C (N)). To see the evolution of the 
pion distribution with increasing N, it is more enlighten- 
ing to study the normalized distribution, defined by 

P(x,N) = H(x,N) J dxxH(x,N), (10) 

where the upper limit of integration is set to 1 on the 
assumption that the contribution from px > Kq is in- 
significant. Thus P(x, N) is the probability distribution 
of producing a tt° at x, for which the differential phase 
space is xdx due to the 2D nature ofpr- 

In Fig. 3 we show P(x, N) for 4 values of N. Note 
how P(x, N) decreases at high x but increases at low x, 
when N is increased. That is the behavior we expect 
when high-pT partons are suppressed, giving rise to low- 
Pt partons. The crossover occurs at around x = 0.06, 
corresponding to px — 0.6 GeV/c. 

Such an evolution of the x-distribution is reminiscent 
of the evolution of the parton distribution in InQ 2 in 
perturbative QCD. Although no precise relationship be- 
tween the two has been established, it is known that in 
the latter case the analytical description is simpler in 
terms of the moments. Thus let us define the moments 

P n (N)= [ dxx n+1 P(x,N). (11) 
Jo 

From Eqs. ©, JjJ) and (JTOJ we can calculate the N de- 
pendencies of P n (N) , which are shown in Fig. 4 for n = 1 
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FIG. 3: Probability distribution P(x,N) for 4 values of N. 



to 5. Evidently, In P n (N) can be well approximated by 
linear dependence on N, i.e., 



lnP n (JV) = a„-b n N. 



(12) 



The slope parameters b n are shown in the inset of the 
same figure. The dependence of b n on n is also linear. 
Thus we may rewrite Eq. I|12jl as 



^P n{ N) = -X 



X = 5.542 x 1(T 4 . (13) 



This is a very economical way of describing the degra- 
dation property of the pion distribution in terms of one 
basic parameter A. 

A physical interpretation can readily be given for A 



when we consider n = 1, for which Pi(N) = (x) N , 
average x at N. From Eq. we obtain 



the 



(x) N = {x) No exp[-A(JV-JV )], (14) 

which exhibits explicitly the exponential decrease of (x) N 
with increasing N, a behavior that solidifies our physical 
notion of what the dense medium does to (pt)- For N — 
2 and N = iV max , we get 



(x)„ I (x) 2 = 0.825, 



(15) 



which gives a quantitative measure of the degree of degra- 
dation. From Eq. (|13fl it is easy also to show that 



il r 

dN 



— yz > N = 



(16) 
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where (x n ) N — P n (N). Hence, the normalized moments 
of P(x, N) are invariant in N. That is a clue to another 
invariant form of the distribution. 

Before we examine the implications of that clue, we 
note that the properties of P n (N) displayed in Fig. 4 and 
described by Eq. (|12fl cannot be expected to be valid for 



FIG. 4: iV dependence of the moments, P n (N), whose log 
values are raised by the quantities in the parentheses. The 
inset shows the slopes b n , the line being a linear fit. 



arbitrarily large n, since the definition of P n (N) in Eq. 
(flip puts more weight on the high end of x when n is 
large. Our cutoff at x = 1, corresponding to px = Kq = 
10 GeV/c, is based partly on the lack of data at higher 
Pt and partly on the recognition that the contribution 
from pt > Kq is unimportant when n is not too large. 
To test the validity of our procedure, we have carried 
out the analysis for Kq = 20 GeV/c, using the same 
&(z), and found that Eq. (|13l) remains to be an excellent 
approximation of the n dependence shown in Fig. 4, and 
that the value of A is larger by just 2%, which is less than 
the experimental errors. Thus we claim that our analysis 
is stable under variations of Kq so long as we consider 
K > 10 GeV/c and n < 5. 

The invariance of the normalized moments in Eq. l|ltj|) 
suggests that we should consider yet another scaling vari- 
able 

u = x/ (x) N = p T / <p T ) N (17) 
for any fixed N. Let us now define 

9(u,N)= (x) 2 N P(x,N), (18) 
whose moments are defined by 



duu n+l ^>(u,N). (19) 



Transforming this integral to an integration over x, we 
find that 

*„(JV) = (x)~ n P n (N). (20) 

It then follows from Eq. lltj|) . that 

d^ n {N) /dN = 0. (21) 

Hence, ^> n (N) is independent of N and we have a scal- 
ing function ^! n , which in turn implies that vE'(m) is also 



4 



independent of N. Indeed, from Eqs. © and (|10|) we see 
that (|18fl can be reexpressed as 

*(u) = $(z(u)) / J duu$(z(u)) , (22) 

where, by virtue of Eqs. JSJ and (|17|) . 

= { X ) N A(N) u. (23) 

Although ( a;)jy A(iV) may appear to depend on TV, it ac- 
tually is a constant 

, . , , f dzz 2 $(z) , , 

7 = (,) JV A(iV) = (,) = L^ = 0.414. (24) 

Thus using z = 7U in Eq. (|22ll . we obtain the scaling 
function \I/(it) 

*(it) = 2.1 x 10 4 (u 2 + 11.65)- 4 ' 8 (1 + 25e~ 1 - 864 ") . (25) 
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FIG. 5: Scaling distribution $(it). 

In Fig. 5 we show ty(u) whose shape differs from that 
of $(z) at low u because of the difference in the power- 
law violating constants, ^(u) is a universal form of all 
the P(x, N) shown in Fig. 3. 



The scaling property of ^>(u) is analogous to the KNO 
scaling of the multiplicity distributions P m (s) in hadronic 
collisions for < 200 GeV 2\. It was found that in 
terms of the scaling variable z = m / (m) , where m is 
the multiplicity, the KNO function ip(z) — (m) P m (s) is 
independent of s. Here, we find that ^(u), defined in Eq. 
(|18|l . is independent of centrality when the scaling vari- 
able, u = pt I (pt), is used. As it is with KNO scaling, 
we have 




1 *(«) = !, for n = 0,l. (26) 



The higher moments are what characterize the scaling 
function, and perhaps scaling violation at some point. 

Recall now that the energy scaling distribution found 
in is the same as the one in Fig. 1 here, although 
the scaling factor K(s) there is different. Since Eq. I|22|) 
is independent of K(s) or K(N), ^>(u) is thus also the 
scaling distribution for any energy. 

It should be emphasized that no theoretical models 
have been used in any part of this investigation. The 
discovery of a scaling behavior over the whole pt range 
that has been measured offers a simple form of the pt 
distribution for dynamical models to describe at any cen- 
trality and energy. The scaling distribution provides us 
with not only a simple picture of the complex px problem, 
but also a way of quantifying the degree of degradation 
of the transverse momentum in the dense medium. More 
importantly, the mere existence of the scaling behavior 
presents a phenomenological obstacle to the realization 
of the theoretical expectation that deconfincment results 
in an anomalous dependence of the pr distribution on 
centrality. The distribution, &(z) or ^>(u), indicates that 
there is no irregularity suggestive of scaling violation, as 
N is varied over the whole range allowed by the Au+Au 
collisions. 
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